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Temporal  and  Spatial  Structure  of  Short 
Wavelength  Infrared  Auroral  Backgrounds 


1.  INTRODUCTION 

During  April  197  9,  September  1979,  and  August  1980,  AFGL/OPR  conducted 
a  series  of  measurements  of  short-wavelength  infrared  and  visible  emissions 
from  the  auroral  particle-excited  ionosphere.  Observations  were  conducted  from 
AFGL's  Optical  Flying  Laboratory  (USAF  NKC-135A  S/N  53120)  during  flights  into 
the  auroral  zone  near  70°  W  longitude  (1979,  from  Pease  AFB,  N.H. )  and  147°  W 
(1980,  Eielson  AFB,  Alaska).  The  principal  goal  of  this  program  is  to  investigate 
the  sources  of  2.  8  —  3. 1  Mm  band  emission  resulting  from  excitation  of  air  by 
energetic -particle  bombardment,  and  to  establish  the  spatial  and  temporal  cor¬ 
relation  of  emission  with  the  deposition  of  ionizing-particle  energy  in  the  upper 
atmosphere. 

The  stimulating  source  of  input  power  is  energetic  (*  5  keV)  electrons  that 
precipitate  into  the  upper  atmosphere  at  high  latitudes  in  association  with  magneto- 
spheric  substorms.  Spatial  distribution  of  flux  and  peak  energy  deposition  altitudes 
of  these  charged  particles  are  measured  from  the  air  fluorescence  they  excite  at 
visible  wavelengths,  that  is,  the  optical  aurora.  In  particular,  fluorescence  in  the 

+  o 

N2  First  Negative  (0-0)  39 14 A  band,  which  is  directly  proportional  to  the  total 
energy  deposited  in  the  atmosphere,  is  monitored  by  a  narrow -field  photometer 
whose  optic  axis  has  been  coaligned  with  a  SWIR  radiometer  having  an  identical 

(Received  for  publication  8  October  1981) 


field  of  view.  These  data  are  supported  by  energy  distributions  measured  by  a 
low  light  level  video  camera,  a  photographic  all-sky  camera  and  a  12 -channel 
filter  photometer  from  which  altitude  prc.'iles  of  energy  deposition  can  be  deter¬ 
mined  (Appendix  II  of  Reference  1)  and  altitude  profiles  of  ionospheric  electron 
density  measured  by  the  Chatanika  radar. 

The  species  ionized  and  excited  by  direct  electron  impact  and  by  secondary 
processes  react  with  ambient  molecules  and  atoms,  producing  chemiluminescence 
that  is  measurable  against  other  natural  IR  atmospheric  emissions.  The  yields  or 
chemiluminous  efficiencies  in  emission  features,  and  the  times  for  buildup  and 
decay  of  the  infrared  radiations,  complement  emission  spectra  in  identifying  the 
radiating  species  and  the  excitation/deexcitation  processes. 

The  atmosphere  above  the  subsonic  jet  aircraft's  operating  altitude  of  35  to 
39  kft  is  essentially  transparent  to  downward-directed  radiation  between  2.  8  and 
3.  1  pm.  A  radiometer  sensitive  to  this  band  and  a  photometer,  each  with  1/6 

radian  fields  of  view  were  pointed  toward  the  zenith  in  flights  through  the  auroral 
2 

oval  in  1975-1978.  In  1979  and  1980  the  fields  of  the  coaligned  instruments  were 
reduced  to  0.  007  radian  (<  1  km  laterally  projected  to  auroral  altitudes)  and  the 
integration  time  required  to  produce  acceptable  signal/noise  ratio  was  lowered 
from  30  sec  to  0.  1  sec  (refer  to  Section  II  for  instrument  specifications).  The 
instrument  configurations  are  compared  in  Figure  1. 

Data  from  the  wide -field  instruments  used  during  the  1975-1978  flights  have 
been  evaluated  and  interpreted  in  References  3  and  4.  The  results,  which  refer 
to  air  between  100  and  150  km  altitude,  can  be  summarized  as  follows: 

a.  SWIR  output  follows  particle -en<  rgy  input  to  within  the  ~20  km,  3/4  min 
resolution  of  the  instruments.  (No  lag  or  "smear"  can  be  resolved  with  confidence 
in  the  data  traces.  ) 

b.  Yield  of  SWIR  photons  is  linear  with  input  at  all  column  input  rates  ob- 

2 

served  (about  5  to  50  erg/cm  -sec). 

c.  Lumped  parameter-efficiency  for  emission  of  what  is  interpreted  as  the 
vibrational  overtone  of  nitric  oxide  molecules  shows  statistically- significant  dif¬ 
ferences  among  both  long  0/2  hr)  and  short  (1  min)  segments  of  data,  with  values 
lying  between  0.  4  percent  and  0.  7  percent. 

1.  Kofsky,  I.  L.  ,  Villanucci,  D.  P.  ,  and  Sluder,  R.B.  (1977)  Data  Reduction  and 

Auroral  Characterizations  for  ICECAP  III.  HAES  Report  No.  59,  DNA 
4220F. 

2.  Huppi,  R.J.  ,  and  Reed,  J.  W.  (1977)  Aircraft  Borne  Measurements  of  Infrared 

Enhancements  During  ICECAP  1975  and  1976,  AFGL-TR -77 -0232, 

AD  A05l  454. 

3.  Kofsky,  I.  L. ,  Villanucci,  D.  P.  ,  and  Andrus,  W.S.  (1978)  Assessment  of 

Rocketborne  and  Airborne  Infrared  Data.  DNA  4731F. 

4.  Kofsky,  I.  L.  ,  and  Villanucci,  D.  P.  (1979)  Further  Assessment  of  Infrared 

Data  from  Aircraft  and  Rocket  Probes.  DNA  5127F. 
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Figure  1.  Comparison  of  Fields  of  View  of  NKC-135  A/55-3120's 
Photometers  and  Radiometers  in  Earlier,  and  in  1979-80  Meas¬ 
urements  of  SW1R  Emission  Structure.  The  diagram  in  the  upper 
right  illustrates  the  aircraft's  trajectory  against  the  rotating  au¬ 
roral  oval  (locus  of  most  intense  particle  precipitation)  in  flight 
on  a  geomagnetic  N-*S-*N  track. 


It  should  be  noted  that  these  aircraft  results  are  by  and  large  consistent  with  those 
from  rocketborne  radiometers,  and  (with  the  exception  of  the  variability)  with 
results  of  laboratory  investigations  of  the  radiative  properties  of  NO. 

These  earlier  flights  also  showed  that  both  SWIR  emission  and  the  referencing 
visible  air  fluorescence  can  be  measured  when  the  auroral  altitudes  are  sunlit  (as 
reviewed  in  Section  5  of  Reference  4).  Navigating  to  compensate  for  winds  and 
operating  near  the  equinoxes,  the  aircraft  is  able  to  hold  the  desired  solar  eleva¬ 
tion  while  flying  within  the  (average)  auroral  oval  and  permitted  airspace.  Veri¬ 
fication  of  a  day-night  difference  in  SWIR  dynamics  and  yield  would  indicate  a  de¬ 
ficiency  in  current  understanding  of  the  processes  that  lead  to  vibrationally  excited 
NO.  No  night-twilight  difference  in  concentrations  of  ambient  reacting  species 


(O,  O 2,  Ng)  has  been  proposed  that  would  have  substantive  effect  on  the  reaction 
chain  that  produces  NO^ .  To  assess  the  behavior  of  SWIR  emission  during  "sunlit" 


aurora,  several  transition-twilight  flights  were  made  with  the  narrow -field  instru¬ 


ments  in  1979  and  1980. 
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The  narrow -angle  instruments  are  tilted  approximately  13°  forward  of  the 
flying  aircraft  to  point  in  the  direction  of  the  geomagnetic  field  when  the  aircraft 
is  within  the  auroral  oval  and  flying  southward.  In  this  mode  they  measure  the 
emissions  from  regions  in  whicr.  the  precipitating  auroral  electrons  have  a  unique 
initial  energy  spectrum.  Accurate  co-alignment  (<  3  arc  min)  is  accomplished 
using  a  precision  optical  bench,  the  design  of  which  is  documented  in  Reference  5. 
A  second,  wider-field,  multi-channel  photometer  (2°  full  angle)  separated  by 
approximately  43  feet  from  (and  nominally  aligned  with)  the  prime  instruments  has 
as  its  major  function  determination  of  energy  deposition  altitudes  from  the  photon 
ratios  of  auroral  features  near  6300A  and  4278A  (Appendix  II  of  Reference  1).  This 
information  is  useful  in  verifying  the  extent  to  which  the  factor  of  2  variability  in 
SWIR  yield  seen  in  earlier  aircraft  measurements  may  be  correlated  with  emission 
altitude.  Calculations  of  collisional  quenching  of  NO1  by  atomic  oxygen  (Refer¬ 
ence  6),  in  fact,  predict  a  minimum  in  SWIR  yield  near  100  km. 

The  determination  of  emission  altitudes  from  measurements  of  the  8300A  and 
4278A  intensities  is  based  on  the  principle  that  this  altitude  is  determined  by  vari¬ 
ations  with  height  of  atmospheric  density  and  composition  that  cause  electrons  of 
different  energies  to  be  stopped  at  different  altitudes.  Major  assumptions  in  the 
theory  are:  (1)  the  optical  axis  of  the  photometer  is  aligned  with  the  magnetic 
field;  (2)  steady  state  excitation  and  loss  of  the  0(*D)  upper  state,  which  radiates 
at  6300A;  and  (3)  a  Maxwellian  spectrum  of  the  energy  of  the  incoming  electrons. 
The  first  of  these  is  generally  satisfied  during  segments  of  flights  when  the  air¬ 
craft  heading  is  near  180°  (magnetic).  In  other  flights,  however,  particularly  the 
twilight  missions  where  the  aircraft  trajectory  was  chosen  to  maintain  a  constant 
solar  depression  angle,  this  condition  is  not  very  well  met.  Thus,  the  confidence 
placed  on  altitudes  derived  by  this  method  depends  (at  least  in  part)  on  aircraft 
heading. 

Additionally,  the  assumed  shape  of  the  particle  energy  spectrum  is  a  potential 
source  of  error.  (The  spectra  measured  in  arcs  from  the  various  ICECAP  rockets 
are  generally  poor  fits  to  the  Maxwellian  model.)  To  this  extent,  the  1980  flight 
series  was  coordinated  with  Chatanika  radar  operations  which  directly  measures 
altitude  profiles  of  ionospheric  electron  densities. 

A  total  of  about  15  hr  of  coordinated  measurements  from  5  data  missions 
originating  at  Eielson  AFB  were  made  in  which  the  aircraft  was  flown  along  (or 
near)  the  magnetic  north-south  meridian  through  Chatanika,  Alaska.  Separate 
flight  profiles  were  calculated  for  each  night  to  obtain  an  adequate  sampling  of 

5,  Sluder,  R.  B. ,  Andrus,  W.S.  ,  and  Kofsky,  I.  L.  (1979)  Aircraft  Program  for 

Target,  Background,  and  Sky  Radiance  Measurements,  A^'IiL-Th -79-0139. 

AD  A 076  959. 

6.  Archer,  D.H.  (1979)  Requirements  for  Improved  Infrared  Prediction  Capa¬ 

bility.  HAES  Report  No.  78,  DNA  4585F. 


both  sunlit  and  nighttime  aurora.  An  attempt  was  made  to  synchronize  the  air¬ 
craft  trajectories  with  the  radar's  scan  rate  so  that  the  13°  tilted  instrument 
fields  would  overlap  the  radar  scan  volume  at  100  km. 

Since  the  radar  facility  requires  a  preprogrammed  tracking  routine  to  attain 
sufficient  precision  for  the  measurements,  no  provision  in  the  radar  track  could 
be  made  for  the  effects  of  winds  on  the  aircraft.  The  aircraft  was  flown  along  a 
track  which  adjusted  for  the  "crabbing"  of  the  aircraft  caused  by  winds  so  that  the 
instruments  would  view  the  same  region  being  viewed  by  the  radar.  The  aircraft 
utilized  an  inertial  navigation  system  (INS)  to  determine  its  position,  but  since  the 
accumulated  error  of  the  INS  from  take-off  to  landing  was  as  much  as  3  km,  the 
exact  overlap  of  radar  field-of-view  and  the  airborne  instrument  field-of-view 
cannot  be  guaranteed. 

It  was  decided  to  set  up  a  series  of  checkpoints  at  which  the  aircraft  would 
arrive  "on  time"  in  conjunction  with  the  radar  motion.  No  course  deviations  were 
taken  from  the  predetermined  track  (about  220  km)  to  "chase"  auroral  forms, 
which  of  course  lessens  the  probability  of  encountering  intense  particle  bombard¬ 
ment.  The  number  of  meridian  cycles  completed  by  the  radar  per  flight  varied 
from  3  l/2  to  7  depending  on  mission  duration. 

There  are  only  two  positions  along  the  track  when  the  aircraft  instruments 
and  the  radar  view  along  the  same  line  of  sight  and  only  one  where  they  are  view¬ 
ing  along  the  same  line  of  sight  and  along  the  magnetic  field  lines,  which  would  be 
the  ideal  condition.  In  the  data  samples  analyzed,  this  ideal  condition  did  not 
occur.  The  altitude  of  peak  energy  deposition  obtained  from  the  radar  measure¬ 
ments  tended  to  be  5  -  15  km  higher  than  altitudes  derived  from  the  aircraft 
12 -channel  photometer  readings  and  were  in  best  agreement  when  the  altitudes 
obtained  by  the  photometric  method  were  averaged  over  a  short  section  of  data 
rather  than  determined  only  at  the  point  of  maximum  auroral  intensity.  This  is 
expected  since  the  radar  data  was  typically  averaged  over  a  15-sec  period.  Radar 
measurements  were  only  available  for  the  1980  flights  from  Alaska,  so  to  standard¬ 
ize  the  results  presented  here,  all  peak  energy  deposition  altitudes  were  derived 
by  the  photometric  method. 


2.  INSTRUMENTATION 

The  measurements  of  visible  and  infrared  aurora  contained  in  this  report 
were  obtained  from  optical  instruments  mounted  in  the  Air  Force  Geophysics 
Laboratory  NKC-135  aircraft.  The  aircraft  provides  an  excellent  platform  for 
measuring  aurora  since  it  provides  not  only  mobility  but  also  the  ability  to  fly  above 
clouds  and  above  the  strong  thermal  infrared  emission  of  the  lower  atmosphere. 
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The  primary  optical  instrumentation  consisted  of  a  SWIR  radiometer,  a  nar¬ 
row  field  photometer,  and  a  low  light  level  television  (LLLTV)  viewing  through 
the  top  of  the  aircraft  13°  forward  of  the  zenith.  This  configuration  enabled  the 
instruments  to  view  approximately  along  the  magnetic  field  lines  when  the  air¬ 
craft  was  in  the  auroral  zone  flying  magnetic  south.  The  SWIR  radiometer  and 
narrow  field  photometer  had  essentially  matched  fields  of  view  (0.  36°  X  0.  36°) 
which  could  be  accurately  aligned  using  an  optical  bench  mounted  onto  the  outer 
surface  of  the  aircraft.  The  optical  bench  and  alignment  technique  has  been  de- 
scribed  by  Sluder  et  al.  The  low  light  level  television  field  was  larger 
(17.  5°  X  23°)  and  contained  alignable  fiducial  marks  to  indicate  the  field  of  view 
of  the  radiometer  and  photometer.  Secondary  optical  instruments  consisted  of 

■7 

a  near  infrared  radiometer  (NIR),  for  monitoring  OH  at  1.7  pm,  a  12-channel 
photometer,  and  an  all  sky  camera.  1  Table  1  gives  a  summary  of  the 
instrumentation. 

2.1  SWIR  Radiometer 

The  SWIR  radiometer  is  a  vacuum  jacketed,  liquid  nitrogen  cooled  radiometer 
with  InSb  detector  and  coated  germanium  and  silicon  refractive  telescope.  The 
radiation  is  modulated  by  a  tuning  fork  chopper  and  filtered  by  a  bandpass  inter¬ 
ference  filter  before  reaching  the  detector.  The  telescope,  chopper,  filter  wheel 
and  detector  are  all  conductively  cooled  by  a  liquid  nitrogen  reservoir  within  the 
vacuum  jacket.  The  telescope  views  through  a  sapphire  window  and  an  optical 
baffle  which  are  cooled  and  kept  frost  free  by  a  continuous  flow  of  cold  nitrogen 
gas  supplied  by  an  auxiliary  liquid  nitrogen  vaporizer.  The  window  and  optical 
baffle  are  exposed  to  the  air  outside  the  aircraft.  Cooling  the  entire  optical  path 
reduces  the  background  flux  on  the  detector  and  significantly  improves  the 
sensitivity. 

The  relative  spectral  response  of  the  radiometer,  as  measured  by  the  AFGL/ 
OPR  Brower  Labs  monochromator,  is  shown  in  Figure  2.  The  absolute  spectral 
responsivity  was  calibrated  using  both  a  bolt-on  coolable  extended  source  and  an 
external  warm  Barnes  12 -inch  by  12 -inch  extended  source.  A  cold  neutral  density 
filter  was  installed  in  the  sensor  during  the  warm  extended  source  test  and  the 
monochromator  test,  which  permitted  the  sensor  to  view  the  ambient  environment 
without  saturating. 


7.  Huppl,  R.J.  (1977)  A  versatile  radiometer  for  infrared  emission  measure¬ 
ments  of  the  atmosphere  and  targets,  Opt.  Eng.  16:465. 
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Table  1.  Instrumentation  Characteristics 


3800  -  6800A  165”  2”  AFT  EK  2745 

EK  2475  Film  FILM 


281  2.85  2.89  293  297  301  305  309  313  317  3.21 
WAVELENGTH (Mm) 


5600 


I  1  1  I  T 

5800  6000  6200 


WAVENUMBER  (cnr1) 

(a)  (b) 

Figure  2.  Spectral  Response  (a)  SWIR  Radiometer  and  (b)  NIR  Radiometer 


2.2  Narrow  Field  Photometer 

The  narrow  field  photometer  simultaneously  measures  3914A  and  5577A  au¬ 
roral  emission  over  identical  fields  of  view.  This  is  accomplished  by  using  a 
common  telescope  and  field  stop  for  both  wavelengths.  After  the  field  stop  the 
radiation  is  divided  by  a  dichroic  mirror.  The  short  and  long  wavelength  radiation 
is  then  filtered  by  separate  interference  filters  before  passing  to  separate  photo¬ 
multipliers.  A  two-axis  adjustment  of  the  field  stop  permits  alignment  of  the 
photometer  field  of  view  with  the  SWIR  radiometer. 

The  spectral  response  of  the  photometer  was  determined  by  illuminating  the 
full  aperture  of  the  photometer  with  collimated  light  from  the  Brower  Labs  mono¬ 
chromator.  The  spectral  response  for  Flights  906-911  and  924-928  (197  9), 
and  Flights  023-029  (1980),  is  given  in  Figure  3.  The  spectral  filters  were 
changed  in  1980  to  provide  a  more  nearly  optimum  spectral  response.  The  ab¬ 
solute  response  for  the  1980  flights  was  determined  by  viewing  an  extended  diffuse 
reflective  surface  painted  with  Eastman  White  Reflectance  Coating  No.  6080  that 
was  illuminated  with  both  an  NBS  calibrated  lamp  (S/N:A299A)  and  a  high  temper¬ 
ature  point  source  blackbody.  Both  techniques  gave  consistent  results. 


(a)  (b) 

Figure  3.  3914/5577  Photometer  Spectral  Response  (a)3914A  Channel,  (b)  5577A 
Channel  —  Flights  906-911,  924-928,  Solid  Line;  Flights  023-029,  Dashed  Line 


Attempts  to  calibrate  the  narrow  field  photometer  in  the  laboratory  for  the 
1979  flights  led  to  erroneous  results.  The  initial  measurement  of  the  spectral 
response  was  found  to  be  incorrect  because  the  full  entrance  aperture  was  not 

O 

illuminated;  a  spectral  leak  was  found  in  the  3914A  filter  which  made  calibration 
with  a  blackbody  source  difficult;  and  results  from  a  collimated  blackbody  source 
and  an  extended  low  brightness  source  did  not  agree.  Because  of  this,  absolute 
calibration  of  the  photometer  for  the  1979  flights  was  determined  by  cross - 
correlating  the  photometer  auroral  measurements  with  simultaneous  auroral 
measurements  obtained  with  the  12 -channel  photometer  which  had  several  years 

Q 

of  consistent  calibration  history. 

The  responsivity  of  the  SWIR  radiometer,  NIR  radiometer  and  narrow  field 
photometer  are  listed  in  Table  2  for  all  the  1979  and  1980  flights.  The  respon¬ 
sivity  for  the  SWIR  and  NIR  represents  the  values  at  the  peak  of  the  spectral  re¬ 
sponse  as  shown  in  Figures  2  and  3  and  have  not  been  corrected  for  any  spectral 
distribution.  The  NIR  radiometer  values  are  for  a  post  amplifier  gain  of  1.  The 
responsivity  for  the  5 577 A  channel  of  the  narrow  field  photometer  is  the 


8.  Kofsky,  I.  L.  ,  Villanucci,  D.  P.  ,  and  Sluder.  R.  B.  (1980)  Evaluation  of 

Infrared  Simulation  Data.  Final  Report  for  Contract  DNA  001-80-C-0010. 
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Table  2.  Instrument  Responsivity 


Flight  No. 

Date 

Narrow  Field  Photometer 

SW1R 

Radiometer 

NIR 

Radiometer 

<Ch2) 

3914 

5577 

906  and  907 

4. 5  X  10‘3 

7. 32  X  10"3 

5.27  X  10’4 

4.  35  X  10'5 

21  Apr  79;  23  Apr  79 

V/kR 

V/kR 

V/kR 

V/kR 

908  -  911 

2. 94  X  10'2 

2.49  X  10'2 

1.  76  X  10*4 

4.  35  X  10‘5 

25  Apr79;  1,2  May79 

V/kR 

V/kR 

V/kR 

V/kR 

924  -  928 

6.  9  X  10'3 

5.  56  X  10'3 

1.  76  X  10*4 

4.  35  X  10'5 

17  Sep  79  -  28  Sep  79 

V/kR 

V/kR 

V/kR 

V/kR 

023 

3.77  X  10‘3 

1.  2  X  10‘2 

1.  76  X  10~4 

3.  76  X  10~5 

7  Aug  80 

V/kR* 

V/kR 

V/kR 

V/kR 

024  -  029 

2.76  X  10'2 

1.  2  X  10'2 

1.  76  X  10'4 

3.  76  X  10‘5 

V/kR 

V/kR 

V/kR 

V/kR 

*This  calibration  is  somewhat  in  doubt. 


responsivity  at  5577 A.  The  responsivity  for  the  3914A  channel  is  corrected  for 
the  integral  of  the  radiation  from  the  3914A  N*  band. 

2.3  Low  Light  Level  Television 

The  LLJLTV  consists  of  an  ISIT  television  tube  and  selectable  lenses  which 
provide  either  a  17.  5°  by  23°  field  of  view  or  a  165°  diameter  all-sky  field  of 
view.  The  spectral  response  is  determined  by  one  of  four  selectable  filters  in 
combination  with  the  S-20  response  of  the  ISIT  tube.  When  the  17. 5°  by  23°  lens 
system  is  used,  an  illuminated  adjustable  reticle  is  used  to  indicate  the  0.  36°  by 
0.  36°  portion  of  the  field  which  is  viewed  by  the  photometer  and  radiometer.  The 
spatial  resolution  with  the  17.  5°  by  23°  optics  is  about  0.  12°  by  0.  12°.  This  sys- 

9 

tern  has  been  described  in  more  detail  by  Sluder  et  al. 

3.  DATA  ANALYSIS 

A  total  of  15  data  segments  from  10  flights  representing  about  7  0  min  of  high 
quality  measurements  were  initially  selected  from  the  data  base  for  analysis. 

The  criteria  used  in  choosing  the  samples  were  as  follows: 


9.  Sluder,  R.  B.  ,  Andrus,  W. ,  Kofsky,  I.  L.  ,  and  Villanucci,  D.  P.  (PhotoMetrics, 
Inc. )  Private  communication. 
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a.  High  signal  to  noise;  generally  IBC  11+  aurora  or  greater 

b.  Fractional  illumination  of  aurora  including  totally  sunlit,  twilight  and 
nighttime  segments  for  investigation  of  SWIR  yield  dependence  on  solar  angle 

c.  Apparent  delay  (or  lack  of  delay)  of  SWIR  emission  with  respect  to 
particle-energy  input,  during  both  day  and  night  auroral  measurements  for  com¬ 
parison  with  solar  angle. 

Calibrated  radiance  traces  from  the  SWIR  radiometer  and  narrow  field  dual 
channel  fluorescence  photometer  for  each  data  segment  are  documented  in  Appen¬ 
dix  A.  Also  included  are  traces  from  the  wide-field  NIR  radiometer  which  moni¬ 
tors  a  region  of  the  OH  Av  =  2  overtone  (1.  66  -  1.  74  Aim)  and  whose  optical  axis 
pointed  in  the  zenith  in  197  9  and  was  nominally  coaligned  with  the  prime  instru¬ 
ments  in  1980.  background  airglow  emission  from  the  hydroxyl  Av  =1  sequence 
which  overlaps  the  SWIR's  bandpass  can  be  scaled  from  the  NIK  readings  and  then 
subtracted  from  the  3  pm  intensities.  The  ratio  of  in-band  OH  fundamental  to 
overtone  emission  has  been  calculated  to  be  about  2. 5:1. 

It  should  be  noted  that  the  narrow  field  photometer  was  a  new  design  whose 
initial  calibration  is  only  indirectly  applicable  to  the  1979  measurements  (see 
Section  2).  Final  calibration  (for  1979)  was  determined  by  cross -correlating  the 
instrument's  raw  voltage  readings  with  absolute  intensities  measured  by  the  air¬ 
craft's  12 -channel  filter  photometer  (Reference  1)  as  described  in  Section  4  of 
Reference  8.  No  instrument  performance  irregularities  were  encountered  in  1980. 

3.1  Chemiluiuinoug  Yields  and  Delay  Compared  with  Solar  Illumination, 

Auroral  Altitude  and  Magnetic  Index 

Table  3  lists  several  auroral  and  measurement  parameters  considered  im¬ 
portant  in  quantifying  infrared  visible  correlations  and  assessing  apparent  chemi- 
luminous  yields.  The  entries  in  Table  3  are  as  follows:  The  first  column  lists 
the  case  number  assigned  to  each  data  segment.  Next  is  flight  data  and  time  of 
meridian  passage  (which  is  needed  to  determine  solar  elevation)  at  the  midpoint 
of  the  data  segment,  whose  duration  is  in  the  third  column.  Aircraft  position, 
fuselage  heading,  and  direction  of  the  flight  trajectory  follows;  the  latter  two 
azimuths  typically  differ  by  a  few  degrees  because  of  crosswinds.  When  the 
fuselage  is  oriented  at  180°  magnetic,  the  coaligned  input-and-output-measuring 
instruments  point  up  the  field  lines. 

The  elevation  angle  of  the  center  of  the  solar  disk  at  the  aircraft's  zenith  and 
(the  directly  applicable  figure)  at  the  latitude  and  longitude  of  the  120  km  altitude 
intercept  of  the  13°  forward  fields  follows  the  solar  declination  (also  included 
because  it  is  an  input  to  the  program  for  calculating  solar  elevation).  Symbols 
D  and  I  indicate  decreasing  (evening  twilight)  or  increasing  (morning  twilight) 
elevation,  to  indicate  whether  the  region  had  been  sunlight  "predosed.  "  The  next 
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Table  3.  Auroral  and  Observation  Parameters  During  SWIR  Enhancements 


column  is  the  hard  earth's  shadow  height  with  no  correction  for  atmospheric  re¬ 
fraction  (see  Figure  27  of  Reference  4).  The  earth's  shadow  height  is  not  an 
accurate  altitude  for  ultraviolet  illumination  because  solar  radiation  is  scattered 
and  absorbed  as  it  passes  tangentially  through  the  earth's  atmosphere.  The  tan¬ 
gent  altitude  at  which  l/2  the  incident  solar  flux  is  transmitted  through  the  (1966 
U.S.  Standard)  atmosphere  is  34  km  at  wavelength  3400A,  45  km  at  3050A,  and 
70  km  at  2380A.  To  a  good  approximation,  the  minimum  auroral  altitude  reached 
by  half  these  ultraviolet  photons  can  be  found  by  adding  "screening"  heights  to  the 
illuminated  height  listed  in  Table  3. 

Column  intensity  ratios  of  optical -auroral  emission  features  at  the  time  in 
the  data  segment  when  the  energy  input  rate  is  a  maximum  are  listed  next.  The 

o  o 

altitude  of  peak  energy  deposition  at  this  time,  calculated  from  the  6300A/4278A 
ratio  and  4278A  intensity  follows.  In  Flights  908  and  909,  before  the  wavelength- 
cycling  photometer’s  field  had  been  directed  13°  forward,  these  altitudes  refer 
to  the  aircraft  zenith;  in  other  flights  they  refer  to  the  excitation  in  the  fields  of 
the  coaligned  0.  36°  instruments. 

The  last  two  columns  list  the  three -hour  local  (K)  and  global  (Kp)  magnetic 
indexes.  Local  K's  refer  to  the  College,  Alaska  sector  and  thus  are  not  directly 
applicable  to  the  1979  flight  series  (conducted  from  Pease  AFB,  N.  H.  ). 

The  data  set  in  Table  3  includes  measurements  on  both  stable  arcs  and  less 
structured  (diffuse)  aurora  that  were  illuminated  to  varying  degrees,  with  direct 
(and  scattered)  solar  radiation.  Fractional  illumination  of  the  emitting  region 
for  each  case  (Table  4)  was  determined  by  first  assigning  a  model  energy  deposi¬ 
tion  altitude  profile  to  each  arc  (or  glow)  based  on  the  peak  altitudes  in  Table  3 
and  following  the  procedure  in  Reference  10.  The  altitude  distribution  of  incom¬ 
ing  particle  flux  was  assumed  to  remain  fixed  for  the  2-3  min  duration  of  the 
measurement  interval.  (For  Cases  15a  and  b,  an  average  peak  altitude  somewhat 
higher  than  that  listed  in  Table  3  was  used  in  the  calculation. )  An  estimate  was 
then  made  for  the  fractional  energy  deposition  rate  lying  above  the  hard  earth 
illuminated  height  defined  in  column  10  of  Table  3. 

Table  4  shows  that  a  sampling  of  day,  night  and  twilight  measurements  are 
included  in  the  selected  data.  Six  cases  were  identified  as  being  totally  sunlit 
(three  each  from  the  1979  and  1980  series).  A  further  six  cases  represent  par¬ 
tially  illuminated  auroras  (again,  three  cases  from  each  flight  series),  while  the 
remaining  three  are  nighttime.  It  should  be  remembered  that  the  entries  in 
Table  4  do  not  reliably  reflect  the  increased  outscattering  of  UV  photons  at  lower 

o 

tangent  heights.  That  is,  a  recalculation  of  Table  4's  fractions  for  2500-3500A 
radiation  would  have  the  effect  of  decreasing  the  fractional  illumination  values. 

10.  Rees,  M.  H. ,  and  Luckey,  D.  (1974)  Auroral  electron  energy  derived  from 
ratio  of  spectroscopic  emissions,  J.  Geophys.  Res.  79:5181. 
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Table  4.  Fractional  Illumination  of  Aurora 
for  the  15  Data  Segments  Identified  by  Case 
Number.  Values  refer  to  the  "Hard  Earth" 
definition  of  solar  illuminated  height  as  de¬ 
fined  in  the  text 


Percent  Sunlit  - 

Hard  Earth  Definition 

Case  Number 

Percent  Sunlit 

1 

0 

6 

0 

Night 

Aurora 

9 

0 

4 

65 

15a 

68 

15b 

45 

3 

53 

14 

15 

5 

47 

11 

100 

7 

100 

13 

100 

12 

100 

8 

100 

2 

100 

To  quantify  spatial  and  temporal  correlation  between  particle  energy  input  and 
infrared  output,  and  to  evaluate  absolute  SWIR  yields,  cross -plots  of  the  limul- 
taneously  measured  3914A  and  3  pm  band  radiances  were  made.  First,  the  data 
traces  were  manually  smoothed,  then  intensities  were  read  at  nominally  5  sec 
intervals.  For  cases  where  the  duration  of  enhancement  was  short  (case  12  for 
example),  the  sampling  interval  was  reduced  to  approximately  1/2  to  1  sec. 

Emission  from  the  hydroxyl  Av  =  1  sequence  which  overlaps  the  SWIR  radio¬ 
meter  bandpass  was  usually  not  subtracted  from  the  3  pm  readings  before  plotting. 
As  the  contribution  from  this  airglow  remains  essentially  constant  during  individual 
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data  segments,  slopes  of  the  cross-plots  (which  are  a  measure  of  the  chemilumi- 
nous  yield)  are  not  expected  to  be  distorted  (or  changed)  by  this  baseline  offset. 

The  projection  of  the  instrument  fields  at  altitudes  where  incoming  particle 
energy  is  absorbed  is  traversed  by  the  aircraft  in  about  4  sec,  and  the  data  smooth¬ 
ing,  as  mentioned  above,  is  of  the  order  of  a  few  seconds.  If  these  times  are  com¬ 
parable  to  (or  longer  than)  the  time  for  aerochemical  reactions  to  excite  SWIR 
radiation,  or  if  production  and  loss  rates  have  reached  equilibrium,  the  cross-plot 
method  is  a  valid  means  for  determining  average  chemiluminous  yields. 

Examples  of  results  from  three  data  segments  representing  day  (Case  12), 
night  (Case  1),  and  twilight  (Case  3)  measurements  on  the  IBC+  aurora  are  shown 
in  Figure  4.  Column  brightnesses  refer  to  1  sec  intervals  (Cases  1  and  3)  on  the 
radiance  traces  which  had  been  previously  filtered  (each  point  on  the  original 
trace  being  an  average  of  11  points  spaced  0.  01  sec  apart).  For  Case  12,  the 
sampling  interval  was  0.4  sec. 

The  all-sky  images  in  the  appendix  (nominal  4  1/2  sec  exposures)  show  that 
the  visible  auroral  distribution  can  best  be  described  as  a  near-uniform  arc  for 
Cases  3  and  12,  while  Case  1  is  a  non -stationary  "breakup".  All  of  the  photographs 
were  printed  with  the  same  exposure  to  qualitatively  indicate  relative  brightnesses. 

Some  of  the  SWIR  enhancements  (Case  12  for  example)  show  obviously  longer 
decay  times  or  peak  offsets  when  compared  to  the  39 14 A  fluorescence  signal  as  is 
evidenced  by  the  loop  formed  by  the  points  in  Figure  4c.  This  effect  is  seen  in 
both  the  twilight  and  nighttime  flights  and  thus  does  not  appear  to  be  directly  cor¬ 
related  with  solar  angle.  (This  point  is  discussed  further  below.  ) 

The  straight  lines  in  Figure  4  are  least-squares  fits  to  the  data  points  (the  R 
values  representing  statistical  correlation  coefficients,  or  "goodness  of  fit"  with 
a  straight  line).  It  has  been  assumed  that  the  3914A  band  fluorescent  intensities 
have  negligible  uncertainty  in  comparison  with  the  3  Aim  band  SWIR  intensities,  and 
that  a  linear  relationship  exists  between  the  two  variables.  The  assumed  linearity 
is  obviated  by  the  curvature  of  points  in  Figure  4c  as  the  fall  time  of  SWIR  radiance 
is  longer  (by  about  5  seconds)  than  that  for  3914A  fluorescence.  However,  the 
seven  circled  points  in  Figure  4c,  which  refer  to  the  first  2  l/2  seconds  of  the 
8  second  measurement  interval,  show  a  high  linear  correlation  when  considered 
separately. 

Slopes  of  the  "best  fit"  lines  can  be  interpreted  as  being  proportional  to  the 

chemiluminous  energy  efficiency  (or  yield)  for  the  nitric  oxide  overtone  sequence 

if  it  is  assumed  that  all  of  the  SWIR  emission  is  produced  by  vibrationally  excited 
2 

NO  populated  by  N(  D)  +  O ^  reactions  (Section  1  of  Reference  3).  In  Figure  4c  the 
fit  to  all  data  points  is  a  reasonable  estimate  to  the  average  NO  yield  over  the 
measurement  interval  as  about  an  equal  number  of  points  lie  above  and  below  the 
line.  (  This  assumption  is  further  supported  by  a  calculation  ratioing  the  "area" 
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Figure  4.  Cross-plots  of  2.8-3  Mm  and  3914A  Intensities  Meas 
ured  by  the  Narrow  -field  Instruments  at  Three  Different  Solar 
Depression  Angles  Representing  (a)  0  Percent  (Case  1), 

(b)  53  Percent  (Case  3),  and  (c)  100  Percent  (Case  12)  Frac¬ 
tional  Illumination.  The  straight  lines  are  uniformly-weighted 
least-squares-fits.  Confidence  limits  refer  to  90  percent 
statistical  certainty  of  the  slopes;  R -values  are  correlation 
coefficients  as  defined  in  Reference  2 


O 

(in  kR-sec)  under  the  SWIR  enhancement  to  that  under  the  3914A  peak  which  shows 
the  average  "integrated"  photon  ratio  to  be  in  agreement  (within  the  limits  of  the 
calculation)  with  the  best  fit  slope.)  The  uncertainties  shown  on  Figure  4  refer 
only  to  90  percent  statistical  confidence  limits  on  the  calculated  slopes  and  thus 
do  not  reflect  other  sources  of  error  such  as  instrument  effects,  window  frost,  etc. 

Following  the  procedure  in  Reference  3,  and  taking  0.42  as  the  fraction  of  the 
NO  Av  =  2  (COCHISE  production  model)  sequence  overlying  the  radiometer  bandpass. 
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we  determined  chemiluminous  vields  from  the  cross -plots  for  each  of  the  IS  seg¬ 
ments  (cases  in  which  the  aforementioned  "area"  method  was  used  are  identified). 
Results,  listed  in  order  of  increasing  solar  depression  angle,  are  shown  in  Table  5. 
The-  column  labeled  PhM  delay  is  a  semi-quantitative  measure  of  the  time  that 
SWIR  output  "lags"  energy  input.  The  values  were  determined  by  visually  com¬ 
paring  peak  offsets  and  apparent  decay  time  shifts  between  3  Mm  band  and  39 14 A 
band  data  traces.  A  more  rigorous,  mathematical  description  of  this  temporal 
correlation  using  Fourier  transforms  is  described  in  another  section,  results  of 
which  are  summarized  under  the  heading  "MB  LAG".  Data  quality  (based  on  signal/ 
noise  of  the  data  streams)  follow;  the  times  in  the  parentheses  refer  to  duration  of 
the  enhancement  interval.  Fuselage  pointing  direction  and  altitude  of  peak  energy 
deposition,  adapted  from  Table  3,  are  in  the  last  two  columns.  The  average  shown 
at  the  far  right  refers  to  peak  emission  altitude  averaged  over  the  observation 
interval. 

The  efficiencies  listed  in  Table  5  show  a  statistically  significant  variability  of 
about  a  factor  of  2  over  both  the  (long  term)  18-month  measurement  period  as  well 
as  among  separate  flights  of  each  series  (short  term),  and  even  within  individual 
flights  (for  example.  Cases  14  and  15  a  and  b.  Flight  027).  A  similar  variation 
was  noted  in  the  data  obtained  by  the  wide-field  instruments  in  1975-1978  (Refer¬ 
ence  3)  where  the  average  efficiency  was  found  to  be  about  l/2  percent.  Here,  the 
mean  yield  appears  closer  to  1  percent,  and  in  some  cases  as  high  as  1.4  percent. 

To  verify,  or  at  least  identify,  probable  sources  of  this  variability,  cross - 
correlation  plots  relating  several  of  the  likely  auroral  and  measurement  param¬ 
eters  to  the  average  overtone  yield  were  made  (Figures  5  and  6).  The  circled 
numbers  refer  to  the  specific  case  numbers  in  Tables  1  and  5.  Figure  5a  shows 
no  apparent  correlation  with  absolute  3914.A  band  level  which  is  not  unexpected  as 
most  of  the  individual  cross-plots  of  3  pm  band  versus  3 9 14 A  band  radiances  are 
generally  linear.  Results  from  the  earlier  wide-field  measurements  support  this 
finding  (Reference  3). 

Figure  5b  and  d  also  indicate  that  SWIR  efficiency  is  independent  of  fractional 
auroral  illumination  and  global  magnetic  index.  Further,  no  correlation  was  found 
between  efficiency  and  aircraft  altitude  (not  shown).  Some  evidence  for  a  depend¬ 
ence  on  the  ratio  of  hydroxyl  fundamental  to  overtone  emission  is  in  Figure  5c, 
however,  the  apparent  trend  [negative  slope;  (-3.2  ±2.2)  of  the  ratio  Av  =  1  photons 
per  Av  =  2  photons  plotted  against  percent  change  in  yield]  is  opposite  to  what  one 
might  reasonably  expect.  That  is,  if  the  OH  ratio  had  been  artificially  increased 
(due  to  the  presence  of  cirrus  clouds  in  the  instrument  field  of  view,  for  example) 
the  result  would  probably  be  higher  efficiencies  (positive  slope)  due  to  increased 
outscattering  of  shorter  wavelength  photons.  It  is  not  clear  whether  frost  on  any 
of  the  viewing  windows  could  produce  the  effects  seen  in  Figure  5c.  In  any  case. 


\ 

ma 


25 


r 


Table  5.  Chemiluminous  Yields  and  Apparent  SWIR  Delays  Derived  From  the  15 
Data  Segments  Listed  in  Order  of  Increasing  Solar  Depression  Angle  and  Compared 
with  Peak  Auroral  Altitude  and  Aircraft  Heading 


Case 

1 

Ave .  Hff. 
(percent ) 

MB 

Ef  f . 
(per¬ 
cent  ) 

f>e  lav 
(sec ) 

MB 

Lag 
(sec ) 

Data 

Qua  1 i t  v 
-S/  N- 

Fuselage 
Point ing 
Di rec t ion 

Aurora  1 

Alt . (km) 
Peak : Ave* 

11 

8.8 

1.36  !  0.13 

2.66 

0 

0.2 

Fair- Good 
(350  sec) 

W( Ferry) 

116 

7 

9.3 

0.70  1  0.09 

0.68 

2-3 

1  .0 

Good 

(  180) 

NW 

160 

13 

9 .5 

0.92  •  0.13 

(0.8  "area” ) 

1  .04 

5 

1  .0 

Good 

(  120) 

S 

119  124 

+/-2 

12 

10.0 

0.57  -  0.10 

{ "area*' ) 

5 

1  .0 

Good 

(8) 

N 

109 

9 

10.4 

1.04  f  0.10 

0.75 

0-1 

1.0 

Fair- Good 
(240) 

S 

130  144 

+  /-3 

2 

10.7 

1.09  0.30 

("area") 

1 .03 

1.5-3 

0.83 

Very  Good 
( 150) 

N 

134 

4 

11.6 

1.36  i  0.17 

1  .22 

2-6? 

0.91 

Good 

(210) 

N 

120 

15 

11.5  Is 

t  half  1 . 1 10  .04 

0.66 

3 

0.91 

Good 

(300) 

s 

120  130 
+/-S 

12.3  2n 

d  half  1 .2  10.40 

1.37 

2? 

0.83 

Good 

(240) 

s 

122  127 
+/-2 

3 

1 1  .9 

1.11  i  0.08 

1  .09 

1-2 

0.83 

Good 

(45) 

N 

i:c 

14 

12.5 

0.51  t  0.08 

0.93 

1-2 

0.91 

Good 

(300) 

N 

1J  122 

V-2 

5 

12.6 

1.09  t  0.08 

1 .03 

3-4 

0.83 

Fair-Good 

(120) 

S 

128  137 

1 

21  .5 

1.37  i  0.08 

1.16 

2-4 

1  .0 

Good 

(  140) 

NE 

120 

6 

23.0 

0,80  1  0.07 

0.55 

0-1 

0.5 

Good 
( 180) 

N 

120 

9 

30.4 

1.30  t  0.21 

1.31 

0-1 

1  .0 

Good 

(80) 

N 

145 

•Altitudes  are  averaged  over  entire  data  segment. 
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AVERAGE  NO  OVERTONE  YIELD,  PERCENT 

Figure  5,  Average  NO  Overtone  Yield  Compared  with  (a)Average  3914.4 
Intensity,  (b)  Fractional  Auroral  Illumination,  (c)  Hydroxyl  Fundamental/ 
Overtone  Ratio,  and  (d)  Global  Magnetic  Index  for  the  15  Data  Segments 
in  Table  3 


there  are  too  few  data  points  in  Figure  5  to  make  a  statistically  compelling  argu¬ 
ment  supporting  the  apparent  anticorrelation.  Additionally,  the  trend  line  is  signi¬ 
ficantly  "weighted"  by  the  two  data  points  from  Cases  12  and  14  which  further  re¬ 
duces  the  confidence  that  can  be  placed  in  the  best  fit  slope. 

Figure  6  shows  that  there  is  a  weak  correlation  between  average  yield  and  the 
altitude  at  which  most  of  the  incoming  electrons  deposit  their  energy.  The  altitudes 
on  the  bottom  half  of  the  figure  are  from  Table  3  and  thus  refer  only  to  that  part  of 
each  data  segment  where  the  3914A  column  brightness  was  a  maximum.  For  the 
five  cases  in  which  the  aircraft  was  heading  magnetic  South,  peak  energy  deposition 
altitudes  were  averaged  over  the  duration  of  the  measurement  interval  and  are 
shown  at  the  top  of  Figure  6. 

The  slope  of  the  regression  line  fit  to  all  data  points  is  3.  6  (or  an  increase  in 
efficiency  of  about  0.28  percent/km)  and  has  a  correlation  coefficient  R,  of  only 
0.  08,  indicating  a  poor  linear  fit.  If  Case  7  is  omitted,  the  slope  becomes  16.  8 
(0.  06  percent/km),  R  =  0.  49,  and  the  points  lie  reasonably  close  to  the  line.  In 
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Figure  6.  Average  NO  Overtone  Yield  for  the 
15  Data  Segments  Compared  with  Altitude  of 
Peak  Energy  Deposition 


both  cases,  the  large  uncertainty  (90  percent  statistical  confidence  limits)  indicates 
that  the  slope  of  the  best  fit  line  could  be  as  high  as  32  (0.  03  percent/km). 

To  further  investigate  this  matter,  "instantaneous"  values  for  efficiency  and 
altitude  were  determined  at  5 -sec  intervals  for  Case  13  during  which  a  uniform 
(sunlit)  arc  was  scanned  by  the  instruments.  The  results,  which  generally  sup¬ 
port  Figure  6,  suggest  an  increase  in  mean  overtone  yield  of  about  0.  15  percent/km 
with  an  uncertainty  of  +0.  18  percent,  -0.  05  percent/km.  This  finding,  however,  is 
in  disagreement  with  the  simple  calculation  shown  on  Figure  7  which  indicates  that 
the  column  integrated  yields  should  be  relatively  insensitive  to  energy  deposition 
altitude. 

As  previously  mentioned,  some  of  the  SWIR  enhancements  have  longer  rise- 
and  fall-times  (up  to  5  seconds)  than  the  corresponding  3914A  traces  while  others 
show  no  apparent  lag.  Determination  of  the  delay  times  from  Table  5  by  visual 
inspection  of  data  traces  is  a  somewhat  subjective  evaluation  procedure  which  de¬ 
pends  in  part  on  the  time  resolution  of  the  data  plots.  As  some  of  the  original 
radiance  traces  were  presented  with  more  compressed  time  scales,  less  certainty- 
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I  3  5  7  9  II  13  15 

INSTANTANEOUS  RATIO  (SWIR  PHOTONS/3914  A  PHOTON) 

Figure  7.  Instantaneous  Ratio  (SWIR/3914A  Photon)  Plotted  Against  Instan¬ 
taneous  Peak  Energy  Deposition  Altitude  for  Flight  026  (Case  13,  1159:30  - 
12  02:53),  and  Compared  with  Simple  Calculation  of  Expected  Altitude 
Dependence 

can  be  placed  on  the  delay  times  assigned  to  those  segments.  However,  since  some 
measurable  offset  of  SWIR  radiance  is  evident  in  almost  all  15  cases,  it  is  unlikely 
that  this  effect  is  an  artifact  of  the  data  presentation.  The  Fourier  transform  tech¬ 
nique  indicates  somewhat  shorter  absolute  lags  with  a  variability  of  a  factor  of  2. 

In  Table  6,  the  offsets  are  separated  into  two  classes  representing  both  short 
(S2  sec)  and  long  (>2  sec)  delays  for  comparison  with  various  auroral  and  meas¬ 
urement  parameters.  The  tabulated  values  refer  to  the  number  of  data  segments 
falling  into  each  category.  It  is  apparent  that  there  is  no  statistically  significant 
"favoring''  of  long  or  short  delays  under  any  heading  as  about  an  equal  number  of 
cases  are  included  in  each  category. 

Of  particular  importance  is  the  lack  of  correlation  with  sun  angle  (or  more 
importantly)  fractional  illumination  of  aurora.  The  cross-plot  in  Figure  8a,  which 
identifies  the  delays  by  the  case  number,  further  quantifies  this  finding.  (Included 
on  the  figure  are  some  of  the  results  from  the  Fourier  treatment  of  the  data.  ) 

These  data  can  be  summarized  as  follows: 


Table  6.  Correlation  of  SWIR  Enhancement  Delay  with  Several 
Auroral  and  Measurement  Parameters 


Figure  8.  Apparent  SWIR  Delay  Compared  with  (a)  Frac¬ 
tional  Auroral  Illumination  and  (b)  Peak  Energy  Deposition 
Altitude  for  the  15  Data  Segments.  Points  marked  <8>  are 
from  the  Fourier  treatment  of  the  data 


—  Cases  100  Percent  Sunlit 

—  Cases  <  100  Percent,  >  0  Percent  Sunlit 

(Range  0.  15  to  0.  68  Sunlit) 

—  Cases  0  Percent  Sunlit 

Note:  "Sunlit"  refers  to  hard  earth  definition  — 
Altitude  in  km  =  0.  98  (SDA  in  degrees)2 


Percent  Sunlit 


PhM  Delay 

£2  sec  >2  sec 


100  Percent 

<100  Percent, 
>  0  Percent 

0  Percent 


2 

3 

2 
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Figure  8b  shows  the  dependence  of  average  SWIR  delay  on  peak  energy  deposi¬ 
tion  altitude.  Again,  no  correlation  between  the  two  variables  is  evident. 

It  should  be  noted  that  some  of  the  5577 A  enhancements  also  appear  to  lag  the 
energy-input  monitoring  3914A  channel  (Cases  1  and  3,  for  example).  However,  in 
all  cases,  the  offsets  are  found  to  be  <  2  seconds  when  the  same  method  of  visual 
inspection  is  applied  to  the  data. 

3.2  3  fjm  and  5577A  Time  Response  and  Hficienrj  Lsing 
Fourier  Transform  Technique 

The  delay  of  5577 A  emission  and  3  /am  emission  relative  to  the  prompt  3914A 
emission  has  been  determined  by  visually  inspecting  intensity  versus  time  plots  of 
the  sampled  and  plotted  data  of  the  three  channels.  A  more  quantitative  method  of 
determining  the  delay,  which  involves  the  frequency  analysis  of  the  three  signals, 
is  presented  here  along  with  a  speculative  chemistry  model  which  is  in  general 
agreement  with  frequency  analysis  of  the  measured  data. 

Short  time  segments  consisting  of  2048  data  points  that  include  the  passage  of 
the  photometer  field  of  view  across  auroral  features  are  digitized.  The  broad 
band  analog  filtering  of  the  three  data  channels  was  carefully  matched  before  the 
samples  were  taken  so  that  no  artificial  relative  delay  was  introduced.  This  dis¬ 
crete  set  is  then  numerically  processed  by  a  Fast  Fourier  Transform  program  on 
a  Digital  PDP-8E  computer  to  give  frequency  transform  spectra  of  1024  complex 
points.  This  is  done  for  all  three  channels,  the  5577A,  3914A,  and  3  pm  pulses 
for  a  common  starting  time.  We  then  ratio  the  transforms  in  the  frequency  domain 
to  determine  the  transfer  functions,  assuming  the  3914A  pulse  reflects  the  time 
dependency  input  to  a  linear  chemistry  system.  The  measured  transfer  functions 
are  then  filtered  and  fit  to  either  a  theoretical  single  or  double  pole  function  by 
selecting  quenching  coefficients  which  give  the  best  fit.  In  general,  we  are  able 
to  derive  a  simple  pole  from  most  pulses,  which  shows  a  surprisingly  similar 
value  around  1.2  sec"*  from  pulse  to  pulse,  for  the  5577A/3914A  transfer  function 
as  well  as  the  3  pm/39 14A  transfer  function. 

Figure  9  shows  a  sample  of  infrared  3  pm  data  measured  on  Flight  909  and 
referred  to  in  this  report  as  Case  3.  This  sample  of  data  has  a  duration  of 
102.40  seconds  starting  at  2  h  13  min  57  sec.  Figure  9  also  shows  the  imaginary 
and  real  components  of  the  Fourier  transform  of  the  data  sample.  The  radian 
frequency  is  proportional  to  the  Kth  transform  point  and  is  given  by 

w  =  2  nK/ (sample  duration  in  sec) 

In  Figure  9,  the  maximum  point  plotted  is  K  =  64  which  corresponds  to  a  frequency 
of  3.  93  radians  per  second.  Figures  10  and  11  are  the  data  samples  and  Fourier 

O  0 

transforms  for  the  3914A  and  5577A  channels. 
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Figure  9.  (Upper  plot)  3  Mm  Auroral  Emission  Measured  2: 13:57  GMT, 
27  April  1979,  Flight  909,  Case  3.  (Center  plot)  Imaginary  compon¬ 
ent  of  Fourier  transform.  (Lower  plot)  Real  component  of  Fourier 
transform.  Time  scale:  12.  8  sec/division,  frequency  scale:  0.  078  Hz/ 
division 
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Figure  10.  3914A  Signal  and  Corresponding  Complex  Fourier 
Transform  of  Data  Measured  2:13:57  GMT,  27  April  1979,  Flight 
909,  Case  3 
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Figure  11.  5  577  A  Signal  and  Corresponding  Complex  Fourier 

Transform  of  Data  Measured  2:13:57  GMT,  27  April  1979, 
Flight  909,  Case  3 


The  ratio  of  the  transforms  of  the  3  nm  divided  by  the  3914A  is  the  transfer 
function  and  is  shown  in  Figure  12.  The  upper  two  plots  are  shown  before  filter¬ 
ing  and  the  lower  plots  are  the  same  data  after  filtering.  The  filtering  consists  of 
a  Fourier  transform-filter-inverse-Fourier  transform  technique  designed  to  pass 
only  the  lower  frequencies,  because  much  of  the  noise  is  more  shot-like  than 
Gaussian.  Figure  13  shows  the  corresponding  unfiltered  and  filtered  transfer 
function  for  the  5 577  A/ 39 14 A  radiation. 

The  filtered  complex  transfer  function  was  used  to  determine  the  best-fit 
quenching  coefficients  and  photon  ratio  of  3  Mm/3914A  and  5577A/3914A.  The  15 
cases  were  analyzed  with  this  technique  and  the  results  are  summarized  in  Table  7. 
Results  for  all  15  cases  are  given  by  Bruce  and  Eastman.  The  delay  of  the  3  nm 

o  o 

or  5577 A  signals  relative  to  the  39 14 A  is  the  inverse  of  the  quenching  coefficients 
and  is  typically  on  the  order  of  0.  8  sec.  The  delays,  as  determined  by  the  Fourier 
analysis,  are  listed  in  Table  5  in  the  column  labeled  MB  Lag  and  are  seen  not  to  be 
in  consistent  agreement  with  delays  obtained  by  visual  inspection  of  the  data  plots. 

The  3  um/3914A  photon  ratios  listed  in  Table  7  assume  that  all  the  emission 
occurs  at  the  peak  of  the  spectral  response  of  the  radiometer.  If  the  emission  is 
assumed  to  be  from  NO  and  the  atmospheric  transmission  and  the  spectral  response 
of  the  radiometer  are  taken  into  account,  then  the  efficiency  can  be  calculated  by 
multiplying  the  3  Um/3914A  photon  ratios  in  Table  7  by  0.  146.  These  efficiencies 
are  listed  in  Table  5  under  MB  EFF.  The  efficiencies  obtained  from  the  Fourier 
analysis  agree,  in  most  cases,  with  those  obtained  from  the  cross-plot  method 
but  disagree  as  much  as  a  factor  of  2  in  some  cases. 

3.2.1  THE  0(*S)  PULSE  PRODUCTION 

A  number  of  chemical  reactions  leading  to  O(^S)  production  have  been  pro¬ 
posed.  They  include  direct  processes  such  as  the  electron  impact  dissociation 
of  molecular  oxygen, 

O0(X3E')  +  e  -  O^S)  +  0(3P)  +  e 
2  g 

and  the  electron  impact  excitation  of  atomic  oxygen, 

0(3P)  +  e  -  0(!S)  +  e  . 

For  the  atmospheric  region  considered  here  these  are  considered  to  be  of  minor 
importance. 


11.  Bruce,  M.  ,  end  Eastman.  C.  (1981)  Frequency  Analysis  of  Visible  and 
Infrared  Auroras,  Utah  State  University  Report  No.  EDL-SRL-81-4. 


Figure  12.  Complex  Transfer  Function  (Unfiltered,  Upper;  Fil¬ 
tered,  Lower)  of  3  ftm  Response  on  27  April  1979.  Data  sample 
measured  2:13:57  GMT,  Flight  909,  Case  3.  Frequency  scale: 

0.  156  Hz/division,  0.98  radians /division 
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*Thla  calibration  Is  somewhat  In  doubt 


The  dissociative  recombination  of  molecular  oxygen. 
o*<x2n  J  +  e  -  0(*S)  +  0(3P) 

2  g 

has  also  been  proposed.  It  requires  the  formation  of  02  first,  and  then  the  follow- 

on  reaction  with  the  electrons.  A  full  description  requires  a  non-linear  chemistry 

12 

involving  ions  of  other  species  as  well.  Recent  estimates  of  its  importance  indi¬ 
cate  it  has  a  minor  role  at  heights  in  excess  of  100  km. 

Finally,  we  consider  the  energy  transfer  from  excited  electronic  states  of 
molecular  oxygen  and  nitrogen.  The  reaction  of  ground  state  atomic  oxygen. 

0(  P),  with  an  appropriate  electronic  state  of  molecular  oxygen,  02,  as 

o*  +  o(3p>  -  o<1s)  +  o2 

13 

has  been  proposed  by  Barth. 

3  +  3 

The  final  reaction,  the  energy  transfer  from  N9(A  L  )  with  0(  P),  as 

»  u 

N„(aV)  +  0(3P)  -  N-OiVl  +  O^S) 

2  U  2  g 

is  probably  the  main  production  source  for  many  auroral  conditions. 

In  addition,  there  is  the  steadier  background  airglow,  probably  due  to  the 

13  14 

Chapman  mechanism,  and  its  variants.  * 

0(3P)  +  0(3P)  +  0(3P)  -  02  (Xj  V  £4)  +  O^S) 

Following  the  deposition  of  energetic  electrons  there  are  many  excited  states 
of  the  major  atmospheric  species  created  by  the  primary  and  secondary  electrons. 
Because  of  its  brightness  and  short  radiative  lifetime,  a  convenient  monitor  of  the 

2  +  p  + 

temporal  behavior  of  the  electron  flux  is  the  N2  (B  —  X  Eg)  first  negative  (0-0) 

band  at  3914A. 

12 

O'Neil  has  found  that  it  takes  20  ion  pairs  to  produce  one  3914A  photon,  so 
we  have  a  relation  between  the  photon  production  rate  and  the  ion  pair  production 
(about  34.  9  ev/ion  pair)  as 

cf'3914(t)  *  (1/2 o)n(t) 

12.  O'Neil,  R.R.,  Lee,  E.T.P.,  and  Huppi,  E.R.  (1979)  Auroral  0(*S)  produL - 

ti^n  and  loss  processes:  ground  based  measurements  of  the  artificial  auroral 
experiment  PRECEDE,  J.  Geophys.  Res.  84:823. 

13.  Barth,  C.A.  (1961)  The  5577  Angstrom  airglow,  Science  134:1426. 

14.  Chapman,  S.  (1931)  Some  phenomena  of  the  upper  atmosphere  (Bakerian 

Lecture)  Proc.  Roy,  Soc..  (London)  A  1 32 : 3 5 3 ■ 
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where is  the  3914A  photon  production  rate  in  photons/sec-cm3,  and  n(t) 
is  the  rate  of  ion  pair  pi  .Auction  (ion  pairs/sec-cm3). 

The  following  set  of  reactions  are  used  to  describe  the  time  dependent  con¬ 
centrations  of  the  species  N(2D),  N2(A).  NO(v),  and  O^S)  where  we  assume  a 


constant  background  of  02 ,  N2,  and  O. 
N^D)  +  02  -  NO(v)  +  O, 

N2(A)  +  0  -  NO(v)  +  N(2D). 

n2ca)  +  o  -n2  +  0(1S), 

NO(v)  +  N  —  N2  +  0(3P), 

O^S)  +  0(3P)  -  2  0(3P). 

o(1s)  +  o2  -  o2  +  0(3P), 
n2(a)  +  °2  -  n2+o2. 

N(2D)  +  O  -  N^S)  +  O. 

N(2D)  +  NO(0)  -  N(4S)  +  NO(0). 
0(XS)  +  NO(0)  -  0(3P)  +  NO(O), 

We  also  have  the  radiative  decay  set, 
0(XS)  -  0(XD)  +  hi/, 

N  (A)  -  N,(X)  +  hi/. 


kj(v)  *  5  X  10" 12  p(v) 

(1) 

1^  -  1.4  X  10"11 

(2) 

kg  =  0.57  X  10"  U 

(3) 

k4  =  2.  5  X  10"11 

(4) 

k5  £  6  X  10"11  e*3o8/T 

(5) 

kg  =  4  X  10"12  e"85°/T 

(6) 

kj  =  4X10'12 

(7) 

k8=4X10"12e-25°/T 

(8) 

k9  =  7  X  10"11 

(9) 

k1Q  =  5  X  10-10 

(10) 

5577A  radiation, 

A  ■  1. 18  sec"1 

A  »  0,53  sec”1 

A^1  (from  Billingsley). 15 


NO(v)  —  NO(v-dv)  +  hi/ 


15.  Billingsley,  F.  P.  (1976)  Calculated  vibration-rotation  intensities  for  NOCX2*) 
J,  Mol.  Spectr.  61:53. 


In  addition,  there  is  the  production  for  N(2D)  and  N2(A)  per  ion  pair,  written 


as 


N(2D) 


2 

is  the  production  rate  for  Nr D)/ ion  pair. 
1  and  1.  5. 


It  has  a  value  between 


€N2(A) 


is  the  production  rate  for  N2(A)/ion  pair. 
0.  5  and  1.  0. 


It  has  a  vaiue  between 


3.2.2  PREDICTING  THE  0(*S)  CHEMISTRY 

From  this  reaction  set,  we  may  write  the  rate  equation  set  as: 


-£-lN(2D)i  =  «  9  n(t)  -  {k.fo,]  +kftEo]  +KQ[NO(0)]}  [N(2d» 
dt  NCD)  1 

+  kglO]  lN2(A)] 

-^■[N2(A)J  =  «N  (A)  n(t)  -  (k7t02)  +  (1^  +  k3)lO)  +  0.53}  [N2(A)] 

^  [0(1S)1  =  k3lO]  [N2(A)1  -  {kg[02J  +k5[Oj  +k10[NO{0)]  +  1.18}  [o(1s>] 


This  is  a  set  of  coupled  linear  time  invariant  first  order  differential  equations, 
driven  by  the  electron  pulse  deposition  function  rj(t).  Such  a  system  is  amenable 
to  solution  by  linear  transform  methods,  especially  the  Fourier  integral  transform. 
The  Fourier  transform  of  the  derivative  is  expressible  in  terms  of  the  Fourier 
transform  of  f(t),  as 

FT  =  iw  F(u) 


Applying  the  transform  to  the  three  rate  equations,  we  have  in  radian  frequency 
space  w, 

^  W2P)  +  k2[°1  tN2(A)1u, 

r)(w)  '  iu  +  a  9  r?(w) 

NCD) 


lO(1S)]ii 

w 

n(w) 


kglo]  [n2(A)1 

iu  +  a  .  r){w) 

OCS) 


where  the  quenching  coefficients  a  are  defined  as 


a  ,  =  {k.lCU  +  kJO]  +  ko[NO(0)]} 

N  rD)  1  2  B  a 


aN2(A)  =  +  (k2  +  k3>  l°l  +  0.  53} 


00(ls)  *  {k6[02>  +  k5l°l  +  k10lNO(0)l  +  18 } 


Thus,  the  transfer  function  for  the  ratio  of  the  557?A  intensity  to  the  3914A 
intensity  is 


A5577  k3t0]  €N2(A) 

Oo)  (i<J+°o<lS))  (iU  +  aN2(A)) 


This  is  a  double  pole  function,  which  has  for  the  w-»0  limit 


<5(0)  Xf*<0> 

5577/  3914 


2°(1.18)«N2(A)k3lOJ 

*  0(lS)°  *2{A) 


We  observe  a  pole  around  1.2  sec  ,  the  radiative  decay  for  the  O^S)  state  with 

no  quenching.  This  is  consistent  with  the  calculation  of  »  ,  from 

OCS) 

a  .  =  (1.  18  +  4  X  10'12  e"850//T  [O  ]  +  5  X  10'10  [NO]) 

ors)  2 

for  most  altitudes  greater  than  100  km,  and  where  we  ignore  oxygen  atom  quench¬ 
ing  in  this  region,  as  recent  laboratory  measurements  suggest.  We  then  predict 
photon  ratio  values  from  1.4  to  4.3  for  reasonable  values  of  [O]  and  [Og] ,  con¬ 
sistent  with  the  range  of  our  field  measurements. 
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3.2.3  THE  3  pm  PULSE  AND  ITS  INTERPRETATION  AS 
CHEMILUMINESCENT  NITRIC  OXIDE 


The  reaction  of  N(^D)  -with  molecular  oxygen  produces  nitric  oxide  with  a  non- 

1  fi 

thermal  vibration  distribution.  This  initial  production  quickly  becomes  redis¬ 
tributed  by  the  natural  process  of  radiative  decay  and  collisional  reactions  with 
other  molecules.  The  complete  radiative  cascade  process  may  be  described  in 
terms  of  a  set  of  rate  equations  of  the  first  order,  written  in  matrix  form  as 

■—  N  =  (A)  (N)  +  P  -  (Q)  <N) 

where  (N)  is  a  column  matrix  for  the  NO  vibrational  levels,  and  (P)  is  the  time 
dependent  production  column  matrix  whose  entries  reflect  the  initial  distribution 
among  vibrational  levels.  The  square  matrices  (A)  (Table  9)  for  radiative 
cascade  and  (Q)  for  collisional  quenching  use  the  Einstein  values  of  NO  given  by 

1C-  19 

Billingsley  and  values  similar  to  those  used  by  O'Neil.  An  expanded  form  of 
the  matrix  set  is  shown  in  Figure  14  with  the  production  used  by  Archer6  in  Table  8. 

We  treat  the  general  matrix  set  by  performing  a  Fourier  transform  with  the 
result 

(iul +  Q  -A)  <N)  =  (P) 

—  ~  ~  u  U 

and  we  assume  that  P  =  n(t)  (p)  where  the  scalar  n(t)  is  the  time  dependence  of 

the  ion  pair  production  used  earlier,  and  p  is  the  fractional  production  of 
1  fi  fi 

Kennealy  as  used  by  Archer.  The  solution  of  the  nitric  oxide  vibrational  levels 
then  becomes  the  inversion  of  a  12  X  12  matrix  acting  on  p  ,  as 

(N)^  =  (iul  +  Q  -  A)'1  (P)  Mu) 

The  predicted  NO  vibration -rotation  band  spectra  for  a  given  Av  sequence  is 
then  determined  by  the  matrix  operation 

£  =  {B(Av)}  <N) 

where  B(Av)  is  the  Einstein  coefficient  matrix  for  a  given  Av. 


16.  Kennealy,  J.  P. ,  DelGreco,  F.  P.  ,  Caledonia,  G.  E.  ,  and  Green,  B.  D.  (1978) 
Nitric  oxide  chemiexcitation  occurring  in  the  reaction  between  metastable 
nitrogen  atoms  and  oxygen  molecules,  J.  Chem.  Phys.  69:1574. 


Table  8.  Einstein  Emission  Coefficients  Avv'  of  NO  and  the  Relative 
Production  for  the  N(^D)  +  02  Reaction  (Billingsley)^ 


Hi 

4V=  1  (  V '  =  V-  1  ) 

<\v=2  (v*  =  v-2 ) 

P(v) 

0 

0.0 

0.0 

0.1044 

1 

10.78 

0.0 

0.1044 

20.43 

0.46 

0.1044 

29.11 

1 .51 

0.1044 

wM 

36.49 

3.10 

0.1044 

5 

42.93 

4.90 

0.1044 

6 

48.54 

7.29 

0.0982 

7 

53.54 

9.63 

0.0982 

8 

57.48 

12.50 

0.0815 

9 

60. S5 

15.68 

0.0418 

10 

62.77 

19.14 

0.0261 

11 

64.15 

22.94 

0.0158 

12 

64.61 

27.15 

0.0120 

E  p(v)  -  1.0 
o 


The  field  measurements  use  a  filter  which  passes  only  a  fractional  amount  of 
each  band,  so  we  must  form  a  row  filter  matrix  F  to  give  the  total  radiation  R(Av) 
through  the  filter  as 

R(Av>  =  (F)<£) 

The  effective  filter  with  atmospheric  transmission  is  designated  "Total  Trans¬ 
mission"  in  Table  10. 

These  operations  were  modeled  using  a  PDP-8E  computer  and  selected  re¬ 
sults  are  shown  in  the  Tables  11  through  14  for  the  zero  frequency  limit,  u  -t  0. 

For  the  6v  *  2  band  at  3  Mm,  there  are  0.2836  photons/  {N^D)  +  02 }  reaction 
available  in  the  zero  quenching  limit,  but  the  filter  used  by  the  AFGL  radiometer 
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a10  a21 

a31 

0  ■(a21+a20) 

a32 

a42 

0  0 

'{a32+a31) 

a43 

0 


"<al2,  ll+a12,  lo'j 


-[Ml 


fq(l)  0  0  . 1 

NO(l) 

0  q(2)  0  . 

NO(2) 

0  0  q(3)  . 

It 

11 

11 

11 

II 

q(12) 

II 

II 

11 

11 

NO(12) 

1)  -q'(2)  0  0 

q'<2)  -q'(3)  0 

0  q'(3)  -q'(4) 


NO(l) 

NO(2) 


q'(12) 


NO<12) 


Figure  14.  The  Rate  Equations  in  Matrix  Form 
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Table  11.  NO  Overtone  {Photons/N(ZD)  +  }  Reaction 


Quenching  Coefficient  a  for  NO(v)  +  M  — ■  NO(O)  +  M 


a{  sec 

-1) 

V  -  v' 

0 

1 .2 

2 

4 

8 

16 

2-0 

0.0167 

0.0145 

0.0132 

0.0107 

0.0074 

0.0042 

3-1 

0.0317 

0.0287 

0.0269 

0.0231 

0.0175 

0.0112 

4-2 

0.0421 

0,0390 

0.0371 

0.0330 

0.0266 

0.0186 

5-3 

0.0444 

0.0419 

0.0403 

0.0368 

0.0310 

0.0232 

6-4 

0.0438 

0.0418 

0.0406 

0.0377 

0.0329 

0.0260 

7-5 

0.0376 

0.0362 

0.0354 

0.0334 

0.0300 

0.0248 

8-6 

0.0285 

0.0277 

0.0271 

0.0259 

0.0237 

0.0202 

9-7 

0.0174 

0.0169 

0.0166 

0.0159 

0.0147 

0.0127 

10-8 

0.0111 

0.0108 

0.0106 

0.0103 

0.0096 

0 .0084 

11-9 

0.0063 

0.0062 

0.0061 

0.0060 

0.0056 

0.0051 

12-10 

0.0035 

0.0035 

0.0034 

0.0034 

0.0032 

0.0030 

Table  12. 

NO  Fundamental  {Photons/N(ZD)  +  02 } 

Quenching  Coefficient  a 

for  NO(v)  +  M  —  NO(0)  +  M 

alsec-1 ) 

V  -  V* 

0 

1.2 

2 

4 

8 

16 

1-0 

0.8788 

0.7004 

0.6085 

0.4422 

0.2608 

0.1222 

2-1 

0.7426 

0.6452 

0.5901 

0.4788 

0.3323 

0.1880 

3-2 

0.6128 

0.5542 

0.5196 

0.4459 

0.3389 

0.2164 

4-3 

0.4958 

0.4595 

0.4375 

0.3890 

0.3141 

0.2191 

5-4 

0.3896 

0.3673 

0.3536 

0.3226 

0.2723 

0.2033 

6-5 

0.2920 

0.2788 

0.2706 

0.2516 

0.2197 

0.1731 

7-6 

0.2092 

0.2017 

0.1970 

0.1860 

0.1670 

0.1378 

8-7 

0.1312 

0.1274 

0.1249 

0.1192 

0.1090 

0.0929 

9-8 

0.0671 

0.0654 

0.0643 

0.0616 

0.0568 

0.0490 

10-9 

0.0364 

0.0356 

0.0350 

0.0338 

0.0316 

0.0278 

11-10 

0.0178 

0.0175 

0.0173 

0.0168 

0.0159 

0.0143 

12-11 

0.0084 

0.0083 

0.0082 

0.0080 

0.0077 

0.0071 

r  " — ~ - 

I 

|  Table  13.  NO  Overtone  (Photons/N(^D)  +  }  Reaction 


Quenching  Coefficient  a  for  NO(v)  +  M  —  NO(v-l)  +  M 


ci(aec~ 

V  -  V1 

0 

1.2 

2 

4 

8 

16 

2-0 

0.0167 

0.0158 

0.0152 

0.0141 

0.0121 

0.0096 

3-1 

0.0317 

0.0306 

0.0299 

0.0282 

0.0254 

0.0212 

4-2 

0.0421 

0.0409 

0.0402 

0.0384 

0.0354 

0.0305 

5-3 

0.0444 

0.0434 

0.0428 

0.0412 

0.0385 

0.0340 

6-4 

0.0438 

0.0430 

0.0424 

0.0411 

0.0388 

0.0348 

7-5 

0.0376 

0.0369 

1.0365 

0.0355 

0.0337 

0.0306 

8-6 

0.0285 

0.0281 

l .0278 

0.0271 

0.0258 

0.0236 

9-7 

0.0174 

0.0171 

0.0170 

0.0166 

0.0159 

0.0146 

10-8 

0.0111 

0.0109 

0..  .OB 

0.0106 

0.0102 

0.0094 

11-9 

0.0063 

0.0063 

0.0062 

0.0061 

0.0059 

0.0055 

12-10 

0.0035 

0.0035 

0.0034 

0.0034 

0.0032 

0.0030 

Table  14. 

NO  Fundamental  fPhotons/Nf^D)  +  C>2 

}  Reaction 

Quenching  Coefficient  a 

for  NO(v)  +  M  —  NO(v 

-1)  +  M 

at  sec" 

-1) 

v  -  v' 

0 

1  .2 

2 

4 

8 

16 

1-0 

0.8788 

0.7916 

0.7425 

0.6429 

0.5070 

0.3566 

2-1 

0.7426 

0.7033 

0.6794 

0.6262 

0.5414 

0.4264 

3-2 

0.6128 

0.5908 

0.5770 

0.5451 

0.4909 

0.4097 

4-3 

0  .4958 

0.4821 

0.4734 

0.4529 

0.4170 

0.3599 

5-4 

0.3896 

0.3808 

0.3752 

0.3618 

0.3377 

0.2980 

6-5 

0.2920 

0.2864 

0.2828 

0.2742 

0.2584 

0.2317 

7-6 

0.2092 

0.2056 

0  .2033 

0.1978 

0.1877 

0.1702 

8-7 

0.1312 

0.1292 

0.1279 

0.1247 

0.1188 

0.1086 

9-8 

0.0671 

0.0662 

0.0656 

0.0641 

0.0614 

0.0566 

10-9 

0.0364 

0.3594 

0.0356 

0.0348 

0.0334 

0.0310 

11-10 

0.0178 

0.0176 

0.0175 

0.0171 

0.0165 

0.0154 

12-11 

0.0084 

0.0083 

0.0082 

0.0080 

0.0077 

0.0071 
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gives  only  0.  1201  photons/  (N(^D)  +  O-  }  reaction  and  0.  0812  photons/  (N(^D)  +  O0 } 

i  - 1  1 
reaction  when  the  NO  quenching  coefficient  is  16  sec  .  The  photon  yield  per 

2  - 1 
N(  D)  +  Og  reaction  as  a  function  of  quenching  coefficient  in  sec  is  shown  in 

Figure  15  for  the  aV  =  2  overtone  band  and  in  Figure  16  for  the  fundamental  band. 

o 

Thus,  the  photon  ratio  for  the  3  urn  to  3914A  radiation  may  be  written  in  the 
zero  quenching  limit  for  100  percent  transmission  through  the  filter  as 


3914 


k,l°2)  E  ,  I 
1  £  nrp)) 


(0. 2836J  _ 

(1720)  ]k1[o2]  +  k8io] 


(5.672)« 


N^D) 


1  +  (0.  8e 


-250/T.  flO 


riQLi 

IM 


For «  o 

Nt  D) 
115  km. 


=  1.  5,  the  ratio  approaches  8.  5  for  altitudes  between  90  km  and 


4.  CONCLUSIONS 

Data  acquired  by  the  narrow -field  radiometer  and  photometer  flown  at  auroral 
latitudes  in  April  1979,  September  1979,  and  August  1980  represent  a  significant 
improvement  in  spatial/ temporal  resolution  of  measurements  on  auroral -particle 
excited  air.  The  low  background  level  and  increased  sensitivity  of  the  liquid  nitro¬ 
gen  cooled  SWLR  radiometer  increases  electrical  bandwidth  by  a  factor  of  300  and 

reduces  the  field  of  view  by  a  factor  of  30  as  compared  to  the  wide-field  instrument 

2 

(with  warm  optics)  used  in  earlier  flight  series. 

Results  from  a  systematic  evaluation  of  15  data  segments,  each  representing 
typically  2-3  minute  observations  on  IBC  II  aurora  of  various  forms,  show  statis¬ 
tically  significant  variability  of  a  factor  of  2  in  "equilibrium"  chemiluminous 
yield.  The  mean  energy  efficiency  of  approximately  1  percent  (if  the  radiation 
is  NO  overtone)  derived  from  the  slopes  of  least-squares  fits  to  the  data,  is  some¬ 
what  higher  than  that  inferred  from  both  the  laboratory  data  (0. 45  percent)  and 

3 

earlier  aircraft  measurements.  Additionally,  some  of  the  SWIR  radiance  en¬ 
hancements  have  comparatively  longer  rise-  and  fall-times  (of  up  to  5  seconds) 

O 

than  the  corresponding  3914A  traces,  while  others  show  no  visually  detectable 
offset. 
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Figure  15.  NO  Overtone  Photons  per  N(2 D)  +  02  Reaction 


Since  the  inferred  chemiluminous  yield  appears  to  change  within  individual 
flights  (as  well  as  among  flight  series)  it  is  unlikely  that  the  variability  is  an 
artifact  of  the  instrument  response  or  calibration. 

Some  evidence  exists  which  suggests  that  "column  integrated"  efficiencies 
may  depend  on  the  altitude  at  which  precipitating  electrons  deposit  their  energy 
(Figures  6  and  7)  as  might  be  expected  from  the  increased  quenching  of  NOJ  by 
O  at  lower  altitudes.  However,  the  trend  inferred  from  the  data  set  indicates  that 
the  mean  yield  is  much  more  sensitive  to  altitude  than  is  reasonably  predicted 
from  first  principles. 

No  statistically  significant  correlation  with  other  auroral  (or  measurement) 
parameters  is  apparent  in  the  data.  Although  preliminary  reduction  of  some  of 
the  narrow -field  data  seemed  to  indicate  a  trend  towards  higher  mean  yields  in 
solar  illuminated  aurora,  the  results  from  the  expanded  data  base  presented  here 
do  not  support  this  finding. 

The  high  spatial  resolution  of  the  radiance  measurements  show  the  SWIR 
column  emission  to  be  closely  correlated  with  air  fluorescence  (to  at  least 
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Figure  16. 


9 

NO  Fundamental  Photons  per  NfD)  +  Oz  Reaction 


5  seconds  and  in  some  cases  less)  both  when  the  aurora  is  pulsating  and  when 
stable  arcs  are  scanned  by  the  aircraft's  forward  motion.  Thus,  the  aerochemical 
reactions  that  excite  2.  8  —  3.  1  pm  band  output  can  be  said  to  "follow"  the  initial 
energy  input  to  within  a  few  km. 
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Appendix  A 

Data  Samples 

This  appendix  contains  15  samples  of  data  which  were  measured  during  flights 
in  1979  and  1980  with  the  SWIR  radiometer,  narrow  field  photometer  and  NIR 
radiometer.  The  date,  time  and  flight  number  of  the  15  samples  are  listed  in 
Table  Al.  Table  Al  also  contains  much  additional  information  which  has  been 
discussed  completely  in  Section  III  of  this  report.  The  table  is  reprinted  here 
for  convenience.  The  case  numbers  correspond  to  case  numbers  referred  to  in 
other  sections  of  this  report. 

Each  sample  of  data  contains  the  time  history  of  radiation  measured  at  3  pm 

O  O 

with  the  SWIR  radiometer,  1.7  pm  with  the  NIR  and  3914A  and  5577A  measured 
with  the  narrow  field  photometer.  The  data  plots  are  followed  by  a  sequence  of 
all-sky  photos  taken  during  the  data  sample.  Relative  to  the  all-sky  photos,  the 
aircraft  is  moving  from  left  to  right  and  the  region  of  the  photo  being  viewed  by  the 
optical  instrumentation  is  slightly  right  of  center. 

In  Cases  11  through  15  all  channels  were  viewing  along  the  same  line  of  sight. 
In  Cases  1  through  9  the  1.  7  pm  channel  (NIR  radiometer)  was  viewing  15°  aft  of 

o  o 

the  3  pm,  3914A  and  5577A  channels.  Case  10  is  missing  since,  for  some  reason, 
the  magnetic  tape  recorder  did  not  function  during  this  time.  In  Case  1  the  cali- 

O 

bration  of  the  3 9 14 A  channel  is  somewhat  in  doubt. 

The  frequency  response  of  the  data  processing  channels  was  kept  as  close  to 
identical  as  possible.  All  channels  were  recorded  on  FM  magnetic  tape  having  a 
frequency  response  of  2500  Hz.  During  data  reduction,  the  3.7  pm  and  1.7  pm 
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Table  Al.  Auroral  and  Observation  Parameters  During  SWIR  Enhancements 


Aircraft  ia  turntn*  to  S  0015:30  -  0017:56 


signals  were  demodulated  with  a  lock-in  amplifier  having  a  12  db/octave  filter  and 
a  time  constant  setting  of  T  e  4  msec.  This  gives  a  complex  frequency  response 
given  by 


A(f)  -  — 12 - 5  (A-l) 

(fo  +  ifr 

where  fo  *  *  39.79  Hz.  The  data  was  then  sampled  at  0.01  second  intervals. 

When  the  data  was  plotted,  a  running  average  of  11  samples  was  taken  to  make  one 
plotted  data  point.  The  running  average  effectively  convolves  the  data  with  a  rec¬ 
tangular  function  which  has  a  width  of  0.  11  sec.  The  resulting  frequency  response 
is  given  by 

B(f)  =  (A -2) 


where  t  =  averaging  time  =  0.  11  sec.  The  total  frequency  response  is  the  pro¬ 
duct  of  A-  1  and  A-2 


H(f) 


fo2  sinjrTf 

(fo  +  if)5  *Tf 


(A -3) 


Substituting  fo  =  39.  79  Hz,  T  =  0.  11  sec  into  A-3,  H(f)  has  a  magnitude  of  0.  707  at 
f  =  3.  97  Hz. 

The  narrow  field  photometer  signals  were  processed  in  an  essentially  identical 
manner.  However,  since  the  photometer  signals  were  DC  coupled  and  did  not  re¬ 
quire  demodulation,  they  were  filtered  by  a  2  pole  filter  which  had  a  rise  time 
matched  to  that  of  the  lock  in  amplifier  used  with  the  3  *im  and  1.  7  nm  signals. 

Equation  A-3  gives  the  frequency  response  of  the  plotted  data  to  temporal 
variations  in  the  aurora.  The  response  of  the  sensors  to  spatial  variations  is 
somewhat  influenced  by  the  temporal  response  but  is  dominated  by  the  field  of  view 
of  the  sensors.  For  a  0.  36°  by  0.  36°  field  of  view  and  assuming  a  range  from  the 
aircraft  to  the  aurora  of  110  km,  the  sensor  would  view  a  section  of  the 
aurora  0.  69  by  0.  69  km.  Assuming  an  aircraft  velocity  of  430  nmi/h  (0.  22  1  km/sec), 
one  field  of  view  is  traversed  in  3.  12  sec.  The  spatial  frequency  response  is  then 
given  by 


G(f) 


(fj  +  if) 


sin  riL  f 
nLjf 


sin  nL.„f 

-ffv2- 


(A-4) 
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where 

f  *  cycles  per  kilometer 
fj  »  179.8  c/km 

*  0.0243  km 

Lg  K  0.  69  km 

G<f)  is  dominated  by  the  last  term  which  has  its  first  zero  at  1.45  c/km. 

The  variations  observed  in  the  data  sample  could  be  either  the  result  of  tem¬ 
poral  variations  or  spatial  variations  which  result  in  a  time  varying  signal  when 
scanned  by  the  motion  of  the  aircraft.  A  detailed  effort  to  separate  the  temporal 
variations  from  the  spatial  variations  has  not  been  performed,  but  the  television 
and  all  sky  television  tend  to  indicate  that  the  variations  are  primarily  temporal. 
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Figure  Al.  Case  1,  Flight  907,  2  3  April  1979 
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Figure  A2.  All-Sky  Photographs  for  Case  1,  Flight  907,  23  April 
197  9.  Total  field  of  view  is  165°.  The  aircraft  motion  is  from  left 
to  right  with  the  instruments  viewing  a  small  region  slightly  right  of 
center.  Each  exposure  is  4.  5  seconds 
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Figure  A7.  Case  4,  Flight  909,  27  April  1979 
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Figure  A 16.  All -Sky  Photographs  for  Case  8,  Flight  926,  25  Septem¬ 
ber  1979 
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Figure  A 18.  All-Sky  Photographs  for  Case  9,  Flight  927,  27  Septem¬ 
ber  1979 
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Figure  A2fi.  All-Sky  Photographs  for  Case  14,  Flight  027,  19  August 
1980 
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Figure  A28.  All-Sky  Photographs  for  Case  15A,  Flight  027,  19  August 
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Figure  A30.  All-Sky  Photographs  for  Case  15B,  Flight  027,  19  August 
1980 
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Attn:  Dr.  H.  C.  Fitz 
P.  Lunn 

Lt  Col  W.  McKechney 

Aerospace  Corporation 
P.  O.  Box  92957 
Los  Angeles,  CA  90009 
Attn:  J.  Crestwell 
P.  Kisliuk 
F.  Simmons 
T.  Stocker 
I.  Spiro 

Electro-Dynamics  Laboratory 
Utah  State  University 
Logan,  UT  84322 
Attn:  D.  Baker 
A.  Steed 

General  Research  Corporation 
Santa  Barbara,  CA  93111 
Attn:  J.  Ise 

Kaman  Nuclear  Tempo 

816  State  Street  (PO  Drawer  QQ) 

Santa  Barbara,  CA  93102 
Attn:  T.  Stephens 

Lockheed  Missiles  and  Space  Co. ,  Inc. 
3251  Hanover  Street 
Palo  Alto,  CA  94304 
Attn:  J.  Kumer 
R.  Sears 
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Mission  Research  Corporation 
P.O.  Drawer  719 
Santa  Barbara,  CA  93102 
Attn:  D.  Archer 

D.  Sappenfield 

PhotoMetrics,  Inc. 

7  Arrow  Way 
Woburn,  MA  01801 
Attn:  G.  Davidson 

I.  Ko.sky 

Physical  Sciences,  Inc. 

30  Commerce  Way 
Woburn,  MA  01801 
Attn:  G.  Caledonia 

R&D  Associates 
P.O.  Box  9695 
Marina  Del  Ray,  CA  90291 
Attn:  F.  Gilmore 
H.  Mitchell 

Space  Measurements  Laboratory 
Utah  State  University 
Logan,  UT  84322 
Attn:  K.  Baker 
D.  Burt 

Space  Data  Corporation 
1333  West  21st  Street 
Tempe,  AZ  85282 
Attn:  S.  Fisher 
D.  Kush 

Stewart  Radiance  Laboratory 
Utah  State  University 
139  The  Great  Road 
Bedford,  MA  01730 
Attn:  R.  Huppi 

J.  Ulwick 

SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 
Attn:  W.  Chestnut 

Visidyne,  Inc. 

5  Corporate  Place 
S.  Bedford  Street 
Burlington,  MA  01803 
Attn:  J.  Carpenter 
T.  Degges 
C.  Humphrey 
W.  Reidy 
H.  Smith 
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